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ABSTRACT 

We study the contribution of thermal and non-thermal processes to the inverse 
Compton emission of the radio galaxy M 87 by modelling its broad-band emission. 
Through this we aim to derive insight into where within the AGN the X-ray, y-ray, 
and VHE emission is produced. 

We have analysed all available INTEGRAL IBIS/ISGRI data on M 87, span¬ 
ning almost 10 years, to set an upper limit to the average hard X-ray flux of 
/(20 — 60keV) ;< 3 x 10“^^ ergcm“^s“^, using several techniques beyond the stan¬ 
dard analysis which are also presented here. We also analysed hard X-ray data from 
Suzaku/VIN taken late November 2006, and we report the first hard X-ray detection 
of M 87 with a flux of /(20 — 60keV) = 10“^^ erg cm“^ s“^. In addition we analyse 
data from Fermi/LAT, INTEGRAL/iENi-N, and Suzaku/yAS. We collected historical 
radio/IR/optical and VHE data and combined them with the X-ray and y-ray data, 
to create broad-band spectral energy distributions for the average low-flux state and 
the flaring state. The resulting spectral energy distributions are modelled by applying 
a single-zone SSG model with a jet angle oiO = 15 °. We also show that modelling the 
core emission of M 87 using a single-zone synchrotron self-Compton model does repre¬ 
sent the SED, suggesting that the core emission is dominated by a BL Lac type AGN 
core. Using SED modelling we also show that the hard X-ray emission detected in 
2006 is likely due to a flare of the jet knot HST-1, rather than being related to the core. 


Key -words: galaxies: active - galaxies: individual: M 87 - X-rays: galaxies - gamma 
rays: galaxies 


1 INTRODUCTION 


Radio galaxies are a subclass of AGN, which display jets 
that are observed with a large angle oi 9 > 10° with re¬ 
spect to the line of sight, enabling a view of both the jet 
and the core. In the unification model of AGN, they are 
thought to be the radio-loud counterparts of Seyfert galax- 
ies ||Antonucci|[i993| |Urry fc Padovai5||1995| [Beckmann &:| 


Shrader 20121. In the case of blazars y-ray emission is ex¬ 


pected to be detected due to the small jet angle and the 
resulting relativistic beaming towards the observer. The rel¬ 
ativistic jet will dominate the emission, and due to Doppler 
boosting this emission can reach into the y-ray and TeV 
range. However, several non-blazar AGN and in particular 
radio galaxies have also been detected in the y-rays (see for 


example the third Fermi/hNT source catalogue; The Fermi 
|LAT GollaborationpOl^, and a few of these sources have 


been observed at energies > 100GeV (Perkins et al.|[2012|. 


The mechanism driving this high-energy emission in radio 
galaxies is still under discussion. 


M 87 is a FR-I radio galaxy (Fanaroff & Riley 1974 
Laing et al.|1983 1, with a central supermassive black hole of 


mass Mbh = (3 — 6) x 10^ Mq (Macchetto et al.]l997 Geb- 


hardt fc Thomas|2009| Batcheldor et al.|2010 1, at a distance 
of 16 Mpc ( |Tonry||1991 1. The jet angle has been estimated 
at 9 = 15° with respect to the line of sight based proper 
motions of the jet features iBiretta et al.|H)^|. 


M 87 has been detected by Fermi/LAT | Abdo et al. 
|2009p , making it the third radio galaxy to be detected 
in y-rays, after Gentaurus A and NGC 1275. M 87 has 
also been detected in the VHE range during flares, e.g. by 
HEGRA ( jAharonian et ar]|2003| ). This source shows vari¬ 
able emission in the soft (< lOkeV) X-ray regime (Harris 
et al. 20091, and extrapolating from high flux states ob¬ 
served by Chandra (/ ~ 2 x 10“'^^ ergcm“^ s“^), a flux of 
/ ~ 10“^^ erg cm“^ s“^ would be expected between 20 and 
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60 keV. Walter & Neronov (20081 reported a detection with 
/ = (8.6 ±1.8) X 10“^^ erg cm between 20-60 keV us¬ 

ing INTEGRAL/lSGKl data. However, we will show that 
this detection can not be conhrmed using the latest soft¬ 
ware and instrument calibration. 

The proximity of M 87 allows us to image the core sepa¬ 
rately from the jet and the diffuse extended emission in sev¬ 
eral wavelengths like radio, optical, and soft X-rays, but not 
in hard X-rays (> 10 keV) nor at higher energies. At energies 
above lOkeV the resolution of commonly used coded-mask 
instruments (e.g. INTEGRAL, Swift/UNT) is of the order of 
10’. Only recently with NuSTAR’s narrow held instrument 
a resolution of ~ 10” can be achieved at hard X-rays using 
grazing incidence optics. 

By studying the spectral energy distribution (SED) it is 
possible to test whether the high-energy emission originates 
from the core, the jet, or from an extended region. Earlier 
SEDs of M 87 have been represented by a synchrotron self- 


Compton (SSC) type model (see for example Abdo et al. 


2009). SSC models are often used to represent blazar spec¬ 


tra, but they have also been applied to other y-ray detected 


radio galaxies, such as Centaurus A (Abdo et al. 2010aI. 


The latter study was inconclusive, though, on the question 
whether the X-ray domain is dominated by the non-thermal 
jet emission, or arises rather from a Seyfert type core (see 
also Beckmann et al.]|2011 ). 

In this paper we present an upper limit on the aver¬ 
age long-term hard X-ray emission of M 87 using 1.7 Ms of 
INTEGRAL IBIS/ISGRI data and different techniques for 
performance enhancement of INTEGRAL IBIS/ISGRI. We 
have also analysed Suzaku data from November 2006, where 
we have detected M 87 for the first time between 20-60 keV. 
This is combined with other data in soft X-rays provided by 
INTEGRAL/JEM-N, 7-rays by Fermi/LAT, historical ra¬ 
dio, infrared and optical emission of the core, and VHE data 
from H.E.S.S. to create an average SED. To understand the 
2006 Suzaku observations we create two simultaneous SEDs: 
one for the M 87 core and one for the bright jet knot HST-1. 
Error values quoted in this paper are at the 1-a level unless 
indicated otherwise. 


2 DATA ANALYSIS 

2.1 INTEGRAL 

In this study we used all available data on M 87 taken by the 
INTEGRAL mission ( Winkler et al.|[2003 l since its launch. 
Observations have been performed in dithering mode with 
pointed observations, so-called science windows, lasting be¬ 
tween 2000 s and 4000 s. The data cover the time from 2003 
to 2011, and most of the observations are taken after 2008. 
We analysed data from the Joint European Monitor in X- 
rays (JEM-X) and from the IBIS/ISGRI imager. JEM-X is 
a coded-mask instrument that consists of two identical co¬ 
aligned telescopes and operates in the 3 to 35 keV band 
( Lund et al.|2003 1. The field of view is circular with a diam¬ 
eter of 4.8 degrees (fully coded) and an angular resolution 
of 3.25 arcminutes. We created images in the energy ranges 
of 3-10 keV and 10-25 keV for each pointing and combined 
these into a single mosaic image for both JEM-X detectors. 
In the 3-10 keV energy band we detect M 87 with a signifi¬ 
cance of 15(7 and a flux of / = 1.6 x 10“^^ ergcm“^ s“^. In 


the 10-25 keV energy band M 87 is not detectable, with a 
3(7 upper limit of / < 1.2 x lO”^'^ ergcm“^ s“^. 

The INTEGRAL Soft Gamma-Ray Imager (ISGRI) is 
part of the Imager on Board INTEGRAL Spacecraft (IBIS), 
which is also a coded-mask instrument. ISGRI is sensitive 
between 15 keV and 1 MeV ( Lebrun et al.|[2C)03 ). The field 
of view is 9° x 9° (fully coded), and the angular resolution, 
limited by the coded mask technology, is 12 arcminutes. 
The total IBIS/ISGRI data set reaches an effective on- 
source exposure time of 1.7 Ms (see Fig. [^. Using the 
standard Offline Scientihc Analysis (OSA) package version 
9 provided by the ISDG ( Gourvoisier et ^|2003 l we cre¬ 
ated images for each pointing in the energy range 20-60 
keV and combined them into a mosaic. At the position of 
M 87 we derived a detection signihcance of 3.8(7 and a flux 
of / ~ 3 X 10“^^ ergcm“^s“^. However, the mosaic image 
shows strong noise features in the vicinity of M 87, and thus 
the detection cannot be deemed as trustworthy. To improve 
the quality of the image we applied several techniques, a 
summary of which can be found in Table 

In order to quantify the quality of the mosaic images, 
we determined histograms of the significance image and the 
root mean square (rms, Srms) of this histogram, where in 
the ideal case the mean should be 0(7 and the Srms = 1. 
Using the firms we can track the improvement of the mosaic 
quality. For the total mosaic created using OSA 9 we derived 
firms = 1.75 (see Tab.[^. To improve the image we started 
by excluding those science windows from the mosaic analysis 
that showed a high noise level, i.e. we removed all science 
windows with firms > 1.2 (~ 4% of the total). The global 
image quality improved, and the M 87 detection is no longer 
signihcant. We have also processed the good science windows 
with the more recent OSA 10 software. We found a detection 
signihcance of 1.32(7 at the position of M 87 and derived an 
upper limit of / < 3.3 x 10“^^ ergcm“^ s“^, within a mosaic 
with a rms of firms = 1.38 (OSA 10 in Table[^. 

We then produced mosaics per revolution (a revolution 
lasts about 3 days and has an effective exposure time of 
~ 200 ks) to evaluate their quality based on the firms value 
of the signihcance map. For the mosaics we set a lower rms 
threshold of 1.1, because the huctuations in the mosaics are 
more averaged out compared to the single science windows. 
Most of the revolutions that have firms <1.1 are within the 
hrst 6.5 years of the mission (rev. < 800, see Fig. [^. This 
is probably due to the evolution of the background (private 
communication). We then combined the images from revolu¬ 
tions where the rms is firms <1.1. The rms of the combined 
mosaic is firms = 1.18 and we derived for M 87 a 3(7 up¬ 
per limit of / < 4.2 X 10"^^ergcm"2s"^ (Tab. OSA 10 
selected revolutions). We also made a mosaic with the com¬ 
bined revolutions with rms firms >1.1, which resulted in an 
upper limit to the flux of 4.9 x 10”^^ erg cm“^ s“^ and the 
mosaic rms of firms = 1.39. 

Lastly we used a technique which required changes to 
the program ii_shadow_ubc, which is part of the OSA 10 
pipeline. This routine performs the background correction 
using the detector images (shadowgrams). To decrease the 
noise in the images we removed pixels from the borders of the 
shadowgrams, since the borders of the detector contain the 
most unpredictable instrumental background. After testing 
several configurations on single revolutions we used a cut 
of the 3 outmost pixels to create 3 mosaics containing all 
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Figure 1. The evolution of the Srms of the significance image of 
the mosaic per revolution. The line shows the cut for Srms = l-l- 
Later revolutions, starting around revolution 800 (2009), have on 
average a higher rms. 


Figure 2. Suzaku/PIN count spectrum between 15 and 70 keV 
with an absorbed power law model fit. The bottom panel shows 
the residuals of the fit in terms of the standard deviation with 
error bars of size one sigma. 


available observations: one mosaic contains the non-modified 
science windows (control), one mosaic where all science win¬ 
dows have been modified, and one mosaic where cuts were 
applied only for those revolutions with a high rms value. 
To quantitatively evaluate these mosaics we calculated the 
rms Srms of the significance map, where we now considered 
only the inner 10° x 10° since the borders of the mosaics 
are unnaturally smooth due to the cutting procedure of the 
shadowgrams, and the area around bright sources in the field 
is excluded. The rms for the mosaic with all science windows 
modified is the lowest, Srms = 1.40, and for the unmodified 
mosaic the highest, Srms = 1.54. However, we found that 
these techniques did not further alter the upper limit on the 
flux of M 87. 

In INTEGRAL observations the systematic er¬ 
rors are dependent on exposure time and on the number of 
sources in the field and their brightness. In the M 87 field the 
number of sources is relatively low and the sources present 
are only moderately bright, therefore limiting the influence 
of systematic errors. In addition, the systematic errors are 
averaged due to the broad energy range and large observa¬ 
tion time. In this field it has been verified that the systematic 
errors are negligible compared to the statistical errors. 


2.2 


Suzaku 


Suzaku has observed M 87 from November 29 to December 2, 
2006 with an elapsed time of 187 ks in HXD nominal point¬ 
ing mode. We analysed data from both the X-ray Imaging 
Spectrometer (XIS, Koyama et al.|20'07 1 and the Hard X-ray 
Detector (HXD, [Takahashi et al.|2007| ). The XIS instrument 
operates between ~0.2-12.0 keV and consists of three sepa¬ 
rate CCD detectors with a field of view of 18’ x 18’ and a spa¬ 
tial resolution of 1.6-2.0 arcmin. The XIS data were repro¬ 
cessed applying the standard event cuts and a spectrum was 
extracted from a circular region with a 60” radius around 
the core. Due to calibration issues the observations between 
0.4 to ~3 keV show strong residuals. The 3-7 keV band dis¬ 
plays thermal signatures such as a tentative iron Ka line at 
6.7 keV (upper limit to the equivalent width of 200 eV), in¬ 
dicating that the emission originates from a hot diffuse gas. 


This is due to the extraction region used that includes both 
the core of M87 and the surrounding hot gas. Since we are 
interested in the non-thermal core emission, we model only 
the 7-10 keV band, where the thermal emission from the hot 
gas is not dominant. An absorbed power-law model yielded 
a fit of Xu = 1-03 (for 150 d.o.f.), with a fixed Galactic hy¬ 
drogen column density Ah = 2 x 10^° cm“^ and a power law 
index of T = 2.3 ± 0.3 (90% error). The extrapolated 2-10 
keV power law flux is / = (2.51 q' 4) ^ 10“^^ ergcm~^ s“^. 

The HXD is a collimated detector with a field of view 
of 34’ X 34’. The detector that consists of two indepen¬ 
dent systems: silicon PIN diodes that operate in the range 
~10-60 keV and GSO scintillation counters that function 
between ~30-600keV. We reprocessed the data applying 
the standard event cuts. Due to the low count rate in the 
GSO band no significant detection could be extracted from 
this detector. For the PIN detector we extracted a signifi¬ 
cant spectrum between 15 and 70 keV (Fig. [^. The spec¬ 
trum can be represented by an absorbed power law with 
a fixed column density Ah = 2 x 10^° cm“^ and a power 
law index of P = 2.8tQ'4 (90% errors) giving a reduced 
xt = 1-17 (for 12 degrees of freedom). The model flux is 
/ = (l.Odl'lg ’j'g) X 10“^^ ergcm“^ s“^ between 20 and 60 keV. 
The combined XIS and HXD spectrum modelled with an ab¬ 
sorbed power law shows a best fit of Xv = 1-03 (150 d.o.f.) 
and a power law index of F = 2.6 ± 0.2 (90% error). While 
the spectral indices of the XIS and HXD differ, they are 
consistent within the 90% confidence intervals. 


2.3 Fermi/LAT 

The Large Area Telescope (LAT, [Atwood et al. 2009 \ aboard 
the Fermi satellite is a pair-conversion instrument that is 
sensitive between 20 MeV and 300 GeV. We have used all 
available data taken between August 2008 and May 2012, 
with a total effective exposure of 30 Ms. We have selected 
source-class-events (P7SOURCE_V6) between 100 MeV and 
100 GeV in a circular region of 30° around M 87. The large 
extraction radius is necessary for the binned analysis we per¬ 
formed. Events with a zenith angle of more than 100° were 
excluded, and we used the standard cuts proposed by the 
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Table 1. Results of IBIS/ISGRI analysis on the field of M 87 applying different data selection criteria. The detection found with the 
OSA 9 analysis is spurious. The count rates have been converted into fluxes assuming a Crab-like spectrum. 


Method 

Mosaic rms Srms 

Detection significance 
M 87 [0-] 

count rate 
[s-l] 

M 87 3cr upper limit 
[10“^^ erg cm“^ s“^] 

Detection significance 
NGC 4388 [0-] 

OSA 9 

1.75 

3.77 

2.34±0.02 

3.0 

148.0 

OSA 10 

1.38 

1.35 

2.20±0.02 

3.3 

113.5 

OSA 10: selected revolutions 

1.18 

1.70 

2.33±0.03 

4.2 

89.4 

OSA 10: borders cut 

1.40 

1.16 

2.20±0.02 

3.2 

111.3 


Table 2. Results of the binned likelihood analysis of the 
Fermi/LAT data in 5 bins between 0.1-100 GeV with a power 
law of photon index F. The errors are given at Icr confidence 
level. 


bin 

energy range 
[GeV] 

flux 

[10~^^ergs cm“^ s“ 

TS 

4 

r 

1 

0.1-0.4 

4.4 ± 1.5 

36 

1.9±0.6 

2 

0.4-1.6 

3.9 ±0.5 

126 

1.9 ± 0.3 

3 

1.6-3.0 

2.1 ± 0.3 

105 

2.7±0.5 

4 

3.0-10 

2.8 ±0.3 

87 

1.47± 0.05 

5 

10-100 

1.6 ±0.4 

12 

2.29 ±0.06 


Fermi team based on the data quality of the events and on 
the instrument configuration. In order to get a reliable result 
for the source of interest, we included in the maximum like¬ 
lihood fitting procedure all sources that had been reported 
in the second Fermi catalogue (2FGL) within 45° around 
the position of M 87. 

Between 100 MeV and 100 GeV we derived a flux for 
M 87 of / = (2.2 ± 0.3) x 10“® phcm“^ s“^ and a power 
law index of 2.16 ± 0.06, with a test statistic TS = 370, 
which corresponds to a detection significance of ~ 19a. Since 
the source is bright in 7-rays, we divided the energy range 
100 MeV to 100 GeV into 5 logarithmic bins to track the 
spectral evolution with energy. The results are summarized 
in Table where we give the significance, flux and power 
law index per energy range. 


3 SPECTRAL ENERGY DISTRIBUTION 
MODELLING 

In order to derive a spectral energy distribution (SED) 
over the whole electromagnetic waveband, we included ra¬ 
dio to UV data from the literature. We used the INTE¬ 
GRAL IBIS/ISGRI upper limit between 20-60 keV, the 
INTEGRAL/IFM-N detection between 3-10 keV and the 
10-25 keV upper limit together with the Fermi/LAT data to 
represent the long-term average emission, where we assume 
that the majority of the low-flux hard X-ray emission is due 
to the core rather than to the jet knots. We combined these 
data with core detections from radio to infrared from NEDQ 
and VHE data from H.E.S.S., which has observed M 87 be¬ 
tween 2003 and 2006 ( Aharonian et al~]|2006 1.To model the 
2006 Suzaku detection we use simultaneous observations, as 
described in Section 3.2. 


http://ned.ipac.caltech.edu/ 


To represent the broad-band SED we applied a one- 
zone synchrotron self-Compton (SSC) model. This model 
assumes an isotropic population of high-energy electrons 
that emit synchrotron radiation followed by inverse Comp¬ 
ton scattering of the synchrotron photons to higher energies 


(Maraschi et al. 19921. In this simplified model the electron 


population is contained in a spherical volume with radius 
R and a randomly orientated magnetic field B. The vol¬ 
ume moves relativistically with a bulk Lorentz factor T(, to¬ 
wards the observer in a jet with an angle 6 to the line of 
sight. The emission is Doppler-shifted with a Doppler factor 
5 = [ri,(l — /lcos6')]“^. The electron energy distribution in 
the jet-frame is assumed to follow a broken power-law, with 
index pi between the minimum energy Emin and the break 
energy Ehv and index p 2 between E\,r and the maximum 
energy Fmaxi 


N{E) = 


kE-P^ if E^in <E <E^, 
kE~^^ if Ehr < E < Fmax 


( 1 ) 


Here, k is the electron normalisation factor, and pi < 3 and 
P 2 > 3. The peak frequencies are dependent on the break 
energy via 


4 2 eB 5 

Vs = - 

3 ' 2'KmeC 1 z 
for the synchrotron peak, where 


'Jbr — 


Ebr 

rrieC^ 


( 2 ) 

(3) 


and for the frequency of the inverse Compton peak (in the 
Thomson regime) 


4 2 

Vie = ^JbrVs 


(4) 


where we assume that the dominant synchrotron power is 
emitted at the peak of the synchrotron branch. Here, 7 is 
the Lorentz factor of the relativistic electrons in the plasma 
blob. 

For some objects, such as bright flat spectrum radio 
quasars (FSRQ), the SSC model does not properly repre¬ 
sent the SED. In addition to the inverse Compton scattering 
of synchrotron photons, external seed photons from e.g. the 
broad-line region, that are Compton upscattered to higher 
energies can contribute to the inverse Compton emission (ex¬ 
ternal Compton component, Dermer & Schlickeiser 19931. 


Otherwise, multiple scatterings should be considered to ex¬ 


plain the Compton dominance of FSRQs (e.g. Georganopou- 
los et al.|2006 1. 


For the modelling of M 87 we have used a single-zone 
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SSC code developed by A. Tramacer^ (to be released soon) 
where the least-square method is used to find the best 
ht of the numerical modelling ([Massaro et al. 2006 Tra- 


macere et al. 2009 20111. This code allows to apply sev¬ 


eral different electron energy distribution shapes, and we 
choose the broken power law shape to compare the re¬ 
sult of the fitting with previous works. In addition to the 
SSC emission there is also a host galaxy component in the 
model, which shows a peak in the optical that is consistent 
with the data. The flux contribution of the host galaxy is 
vfhost = (4 ± 2) X 10“^^ ergcm“^ s“^ for all hts presented 
below. In the following sections the results of the modelling 
is described (see Tablefor the best-fit results). 


3.1 SSC model for the average low-flux state 


At first, we kept the angle 9 and the bulk Lorentz factor 
Ti, free with a fixed beaming factor S — 5. This value is 
consiste nt with the apparent m otion of ~ 0.5c observed in 
the jet ( Kellgrmann et ah [2007 1. The fit has a of 2.8 (17 
d.o.f., Table|^. This gives a value of Smin = 2.6 x 10® eV for 
the minimum energy of the electrons, JSmax = 6.5 x 10^® eV 
for the maximum energy and Ehr ~ 1.3x 10® eV for the break 
energy. The magnetic field has a value of 13 = 1.7 x 10~® G 
and the radius of the emission region R = 3.3 x lO'^® cm. The 
indices of the broken power law are pi = 1.1 and P 2 = 3.5. 

Next we have hxed the jet angle at 0 = 10°, a value that 
is closer to that of blazars, although the apparent angle as 
seen in the large scale radio jet appears to be rather 9 = 15°. 
This allows us to compare our results of the SED fitting with 


the previous work of Abdo et al. (20091, who also used 9 — 


10° for their SED fitting, although this value does not use 


the true orientation of the jet (e.g. Meyer el^l. 20131. The ht 
gives a reduced x? = 3.2 (16 d.o.f.. Table [^. The magnetic 
held is B = 2.2 x 10~® G and the radius of the emission 
region R — 2.4 x 10^® cm. A bulk Lorentz factor Tj, = 3.4 is 
found, which is slightly higher than the bulk factor derived 


by Abdo et al. (20091 of T;, = 2.3. The minimum energy 


is Emin = 2.1 X 10 eV, the maximum energy is Bmax = 
5.1 X lO'^^eV, the break energy Ehr = 1.3 x 10® eV, and 
the power law indices found are pi = 1.1 and p 2 — 3.5, all 
consistent with the previous result using a free angle. Except 
for the minimum energy of the electrons these values are also 


consistent with the values used by Abdo et al. (20091. 


After this we htted the data with a more conser¬ 
vative angle of 0 = 15°, consistent with proper motion 


observations by Biretta et al. (19991 and Meyer et al. 
(20131. The SED is presented in F^. The fit, with a re¬ 
duced Xv = 1-9 (17 d.o.f.. Table Im, yields a lower mag¬ 
netic field, B = 2.0 x 10“® G, and a larger emitting region, 
B = 5.6 X 10^^ cm, compared to the previous fits. While the 
second index of the broken power law is similar with p 2 = 
3.6, the first index changes to pi = —1.8. The bulk Lorentz 
factor increased slightly to T;, = 3.8. The minimum energy 
of the electrons decreased to Emin — 2.8 x 10^ eV and the 
maximum energy decreased slightly to Emnx = 2.3 x 10^® eV. 
Also the break energy of the electron energy distribution de¬ 
creased to Ehr = 2.0 X 10® eV. 


^ http://www.isdc.unige.ch/sedtool/ 


e 



loglK) (Hz) 


Figure 3. SED of the average low-flux state of M 87, using a 
jet angle of 9 = 15°. The arrows show the INTEGRAL upper 
limits. The solid line from the radio to VHE domain shows the 
SSC fit. The synchrotron contribution, dashed, is dominant at low 
frequencies and the inverse Compton contribution, also dashed, 
is dominant at high frequencies. The host galaxy component is 
shown with a solid line in the range 10^^ — 10^® Hz. The Chan¬ 
dra/ NCf 8 observation of the M 87 core at 2.4 X 10®^ Hz was taken 
at the end of 2008, when M 87 was in a quiescent state (|Ahdoj 
|et al.|2009| |. 


3.2 SSC model for the high flux state in 2006 


Since the Suzaku/RW data taken in 2006 indicate a flux 
level that is 3 times as high as the upper limit determined 
from the average 2003-2011 data provided by INTEGRAL 
IBIS/ISGRI, and similarly the Suzaku/yAS data show a 3-10 
keV flux about four times that of the INTEGRAL JEM-X 
core flux, we also collected additional 2006 data in order to 
derive a simultaneous SED for this period. 

Because we are not able to resolve the core and jet with 
Suzaku/RYN, it is not clear whether the flux increase origi¬ 
nates in the core or in one of the jet knots (see e.g. the light 
curves presented by | Abramowski et al.|20T^ . Therefore we 
create two different SEDs: one for the M 87 core, and one for 
the bright jet knot HST-1, which is known to flare. We add 
simultaneous observations for both components from HST 
( Perlman et al.||2011 (, VLA ( Harris et al.|2069 l and VLB A 
( Cheung et al.|2007 1. In addition we add simultaneous VHE 
data from MAGIG ( |Berger||2011[ ), however also in gamma- 
rays it is not possible to distinguish between the core and 
the HST-1 knot. Lastly, we also include the HESS observa¬ 
tions from the average SED ( [Aharonian et aL]|2006[ ). While 
these observations were not taken simultaneously, in 2006 
no increased emission was observed from M 87. 

The SED for the M 87 core is presented in Fig. The 
fit has a reduced Xv = 3.1 (2 d.o.f.. Table |^, where the 
maximum energy has been fixed to Bmax = 2.5 x 10®® eV. 
Using an angle of 9 = 15°, the bulk Lorentz factor is consis¬ 
tent compared to the low-flux average state with F;, = 3.6. 
The radius of the emitting region of i? = 5.0 x 10®®" cm is 
consistent with the low-flux state modelled with the same 
jet angle, whereas the magnetic field slightly increased to 
B = 3.3 X 10 ® G. Both the minimum energy and the break 


































































6 de Jong et al. 


Table 3. The results of the SED fitting of M 87. The first three columns show the results of fitting the average low-flux state using 
several different configurations for the jet angle. The columns ”2006 core” and ”2006 HST-1 knot” refer to the 2006 SEDs of the M 87 
core and HST-1 jet knot, respectively. The last column shows the result of the SED model presented by |Abdo et al.| ( |2009| |. 


Parameter 

9 = 15° 

Fixed beam 

O 

o 

II 

2006 core 

2006 HST-1 knot 

Abdo 2009° 

e 

15°“ 

_b 

10 °“ 

15 ° a 

15°“ 

10 °“ 

B[G] 

2.0 X 10“® 

1.7 X 10“® 

2.2 X 10“® 

3.3 X 10“® 

6.4 X 10“® 

5.5 X 10“® 

R[cm] 

5.6 X 10®'^ 

3.3 X 10®® 

2.4 X 10®® 

5.0 X 10®®’ 

1.0 X 10®® 

1.4 X 10®® 

Lb 

3.8 

_b 

3.4 

3.9 

1.2 

2.3 

EmiJeV] 

2.8 X 10^ 

2.6 X 10® 

2.1 X 10® 

6.7 X 10'® 

6.0 X 10'® 

5 X 10® 

EmaxjeF] 

2.3 X 10^® 

6.5 X 10®® 

5.1 X 10®® 

2.5 X 10®® “ 

2.5 X 10®® “ 

5 X 10®® 

EbAeV] 

2.0 X 10® 

1.3 X 10® 

1.3 X 10® 

5.0 X 10® 

1.0 X 10® 

2 X 10® 

Pi 

-1.8 

1.1 

1.1 

2.8 

2.5 

1.6 

P2 

3.4 

3.5 

3.5 

3.6 

3.4 

3.6 

xl 

1.9 (17 d.o.f.) 

2.8 (17 d.o.f.) 

3.2 (16 d.o.f.) 

3.1 (2 d.o.f.) 

1.9 (2 d.o.f.) 

- 


a) parameter fixed; 

b) beaming factor 6 fixed to 5, see text; 
c) see |Abdo et al.| | |2009[ l 



KaM (Hz) 

Figure 4. SED of the core of M 87 during the high-flux state in 
2006, using a jet angle of 0 = 15°. The SSC fit is shown with a 
solid line. The synchrotron branch is shown with a dashed line at 
low frequencies, and the inverse Compton contribution is shown 
with a dashed line at high frequencies. 

The Suzaku/PIN data do not seem to match the SED. 

energy of the electron energy distribution increased slightly 
to Smin = 6.7 X 10^ eV, and -Ebr = 5.0 x 10* eV. The indices 
of the power law describing the electron energy distribution 
have indices pi = 2.8 and P 2 = 3.6. As can be seen in Fig. 
the model is not consistent with the Suzaku spectrum, which 
dictates a steep power law in the X-ray regime. Therefore, 
the Suzaku spectrum implies that the hard X-ray data de¬ 
scribe rather the tail of the synchrotron branch than the 
inverse Compton branch. 

Since the jet-knot HST-1 is known to flare we also model 
it as a potential origin of the hard X-ray emission detected 
by Suzaku, see Fig.^ A jet angle of 0 = 15° is used and the 
ht has a reduced Xu = 1-9 (2 d.o.f.. Table [^. Also in this 
case the maximum energy is fixed to i?max = 2.5 x 10^* eV 
to allow convergence. The bulk Lorentz factor is lower com¬ 
pared to the 2006 core fit and the average low-flux state with 
Fi, = 1.2. The size of the emitting region is also smaller, with 



log{t/) (Hz) 

1.0| -^^^^^^-----,- ^^^^^-r 

0.5 - + * - 

^ 0.0---_- 

^ -0.5 

-1.0 ^ 

-1.5 I- ^^^^-'-'-^^^-'-'-'-^ 

9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 

logfw) (Hz) 

Figure 5. SED of the HST-1 knot of M 87 during the high- 
flux state in 2006, using a jet angle oi 6 = 15°. The SSC fit is 
shown with the solid line. The synchrotron branch, dominant at 
low frequencies, and the inverse Compton contribution, which is 
dominant at high frequencies are shown with dashed lines. 

The Suzaku/PIN data match the SED. 

R= 1.0 X 10^® cm. The magnetic field strength increased to 
B = 6.4 X 10“^ G. The minimum energy of the electrons is 
similar to the fit of the core with i?min = 6.0 x 10^ eV, and 
the break energy decreased to Ebi = 1.0 x 10® eV. The power 
law indices that characterise the electron energy distribution 
are similar to the 2006 core ht with pi = 2.5 and p 2 = 3.4. 

As can be seen in Fig. the Suzaku data points lie on the 
inverse Compton branch, rather than on the synchrotron 
branch, consistent with the hard X-ray spectrum. 


4 DISCUSSION 

In the following we first discuss the possible origin site of the 
emission seen in INTEGRAL/JEM-yi., before we turn to the 
hard X-ray variability and the spectral energy distribution. 
Finally, we discuss the possible source type at the origin of 
the high-energy emission in M 87. 
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4.1 JEM-X: core emission 

Due to the the angular resolution of JEM-X (FWHM = 3 
arcmin), the flux of M 87 observed by JEM-X is the sum 
of core, jet and extended diffuse emission. Here we are in¬ 
terested only in the core emission, because we want to in¬ 
vestigate the origins of the high-energy emission, which is 
emitted from, or close to the core. In addition, the mea¬ 
sured JEM-X flux is an average over several years. The jet 
knot HST-1 is known to be variable ( Harris et al.|20d9 l, and 
was in outburst between 2003 to 2007, with the luminosity 
of the knot peaking in 2005. For a light curve of the X-ray 


emission of HST-1 see Fig. 1 in Abramowski et al. (20121. 


Since the JEM-X observations we used are taken mostly af¬ 
ter 2008, and none during 2005, we assume that we did not 
include data during which the knot was in outburst. 

Using literature values from observations of the knot 
HST-1 made with Chandra/ACIS by Perlman & Wil- 


(20051 and by [Wilson fc Yang (20021 for other 


jet knots, we find that the combined jet flux is about 
/(3 — lOkeV) = 10~^^ ergcm“^ s“^. In quiescence the jet 
emission is dominated by the jet knots HST-1, A, and D. 
For the extended emission we used observations from Chan¬ 
dra/ ACIS, taken beginning of May 2005. We analysed the 
data, with a total exposure of 123 ks, and extracted a spec¬ 
trum in a circular region with a radius of 2.5 arcminutes 
centred around the core of M 87. The core and jet were 
excluded in the analysis. This spectrum yielded a flux of 

This implies that 
with 


/(3 — lOkeV) = 9 X 10 erg cm ^s 


the core emission is about / = 6 x 10 erg cm ^ s ^ 


a variability by a factor of two in X-rays (Hilburn & Liang 


2012 ). 


The spectrum extracted from the 2006 Suzaku/'KIS ob¬ 
servation also showed a thermal and non-thermal component 
due to the extraction radius covering both the core and sur¬ 
rounding medium. Since a spectrum is available, disentan¬ 
gling these components is a more straightforward task and 
we find a 3-10 keV flux of / = 2 x 10“^^ erg cm“^ s“^ for 
the non-thermal component, which is almost four times the 
JEM-X core flrrx. The strongly increased flrrx derived from 
the XIS data implies a flare when compared to the average 
flux. 


4.2 Hard X-ray variability 

We found an overall 3 (t upper limit to the M 87 flrrx 
of /(20 — 60keV) < 3 x 10“^^ ergcm”^ s“^ using all 
available INTEGRAL data for a total effective expo¬ 
sure time of 1.7 Ms. Apart from a likely spurious de¬ 
tection by INTEGRAL at 5.1 (j, corresponding to a flux 

■s 


of /(20 — 60keV) = (9 ± 2) x 10 erg cm ^ s ^ ([Walter & 


Neronov] 20081, M 87 had not previously been detected in 


hard X-rays. Using the same data set as [Walter fc Neronovj 
( [2008[ ) but applying the latest analysis software OSA 10 
we derived a Scr upper limit flux of /(20 — 60keV) < 
8 X 10“^^ ergcm“^ s“^, inconsistent with the earlier detec¬ 
tion claim. 

This does not mean, however, that the hard X- 
ray flux of M 87 never exceeds this average value. 
We found that, during a Suzaku/RYN observation at 
the end of November 2006, M 87 is detected with a 
flux of /(20 — 60keV) = 10“'^^ ergcm“^ s“^. There is no 


INTEGRAL IBIS/ISGRI observation at this time to con¬ 
firm the detection, but in July 2006 INTEGRAL has ob¬ 
served M 87 for 23 ks constraining a 30- upper limit of 
/(20 — 60keV) < 3 x 10“^^ ergcm“^s“^, consistent with 
the flux measured by Suzaku/RW. 


4.3 Spectral energy distribution 

To characterise the broad-band emission of M 87 we have 
produced three SEDs; one representing the average low-flux 
state and two representing the increased X-ray emission as 
observed by Suzaku at the end of 2006. The data used in the 
average low-flux SED have been taken at different times, 
which might result in a biased measurement in case of vari¬ 
ability since several spectral states might contribute to the 
data. 

To investigate whether there is contamination of differ¬ 
ent spectral states in the M 87 SED, we consult light curves 
of M 87 from 2001 to 2011 based on radio, optical. X-rays 
and VHE observations ( Abramowski et al.|2012 |. The VHE 
band shows increased activity in 2005, 2008 and 2010, dur¬ 
ing which the X-ray emission also increased by a factor of 3. 
During the VHE flare in 2008 the radio emission from the 
jet base increased as well, indicating the VHE emission is 
likely produced near the SMBH. However, during the VHE 
flare in 2008, no enhanced radio emission has been observed. 
The VHE observation used in the average low-flux SED was 
taken in 2004 when M 87 was in a low state, so there is no 
contamination from any of the observed VHE flares. The 
radio and optical detections of the core used in the SED are 
averaged over several years, and the light curves show varia¬ 
tions by a factor of ~ 2. Variability at this level will not sig¬ 
nificantly alter the physical parameters derived from mod¬ 
elling the SED. A similar approach of using time-averaged 


SEDs has been applied already successfully to M 87 (Abdo 
et al.[[2009 ) and other radio galaxies, for example Pictor A 
( [Brown fc Adams|2012p and 3C 111 ( [de Jong et al.|2012p . 

The high flux detected by Suzaku/RTN indicates a flar¬ 
ing state of the source. The phenomenology of the flaring ac¬ 
tivity of M 87 hints at a behaviour similar to that observed 
in blazars, and in particular in BL Lac objects. During a flare 
an increase of the peak energies of the SED correlates with 
an increase in the corresponding peak fluxes (see for exam¬ 
ple Mrk 501, Tavecchio et al.|2001 1. In 2005, the increase in 
VHE, X-ray and optical emission coincided with the expul¬ 
sion of the jet knot HST-1 from the core. During the Suzaku 
observations at the end of 2006 the light curves presented in 


Abramowski et al. (2012) do not show increased activity in 


the VHE or optical band. Even though the amount of simul¬ 
taneous observations is sparse, based on the 2006 SED mod¬ 
elling the HST-1 knot is the more likely candidate for the 
hard X-ray emission. The SED models show that the hard 
X-ray spectrum can not be reconciled with the simultaneous 
core observations. If the hard X-ray band would be on the 
tail of the synchrotron branch for the 2006 core observations, 
this would indicate a synchrotron peak frequency shift from 
~ 10^^ Hz for the low-flux average state to > 10^® Hz for the 
flaring state. A peak frequency shift this large has not been 
observed even in the most extreme blazar flares. Chandra 
had also observed M 87 on 29 November 2005, but due to 
pile-up no useful spectra could be extracted during this time 
(Harris et al. [20091. Since Chandra/ACIS has a small angu- 
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lar resolution of 0.5”, it is possible to monitor the core and 
jet knot separately. Because of the pile-up during the ob¬ 
servation, the intensity is expressed in detector-based units 
of keV s“^ where the entire energy band from 0.2-17 keV is 
integrated. For the core an intensity of 0.51 ± 0.02 keV s“^ 
was measured, whereas the HST-1 knot has an intensity of 
4.03 ± 0.04 keV s“^. This supports our conclusion that the 
X-ray flare detected by Suzaku originates in the HST-1 knot 
rather than the core. 

Since similar SED modelling has been applied to other 
7 -ray bright radio galaxies, comparing the physical parame¬ 
ters derived from these models with the average fit of M 87 
will help us put the results in a broader frame-work. Due 
to the degeneracy of several SED parameters, we can qual¬ 
itatively compare the bolometric luminosity with the mag¬ 
netic field B, the radius of the emitting region R and the 
Doppler factor 5, as these parameters influence the overall 
SED power. 

The FR-I radio galaxy Cen A has also been observed in 
the gamma-ray band, and similar to M 87 the overall SED 


can also be modelled by a simple SSC process ( Abdo et al. 
|2010b[), although a strong Seyfert contribution is visible in 


the X-ray and optical domain |Beckmann et al.|2012 1. Abdo 


et al. (2010b I present several SED fits to the core of Cen A. 


Comparing the results for M 87 (jet angle 0 = 15°) to the 
model for Cen A (jet angle 0 = 30 °), M 87 displays a slightly 
higher Doppler factor of <5 = 3.9, compared to the Doppler 
factor used for Cen A (5 = 1) due to the smaller jet angle 
6. A lower Doppler factor causes the emission to appear less 
boosted. The magnetic field B and radius R of the emitting 
region are quite different for both sources: for Cen A a value 
of B = 6.2 G was found for the magnetic field, and a much 
lower value of B = 2.0 x 10“® G for M 87. Since the syn¬ 
chrotron emission depends on B, a stronger magnetic field 
will result in a higher synchrotron flux. The radius of the 
emitting region is reported to be R — 5.6 x 10^^ cm for M 87 
and i? = 3 X 10^® cm for Cen A (Abdo et al. 2010b I. The 


radius R defines the amount of radiating particles, therefore 
a larger emitting region results in a higher flux. The lower 
magnetic field strength B used to model M 87, combined 
with the larger radius R compared to Cen A then results in 
a similar overall power in the SED, consistent with the com¬ 
parable bolometric luminosity of about Lboi ~ lO^^ergs”^ 
of the two sources. 

Comparing the SED parameters with the luminous 
(Bboi = 5 X lO'^'^ergs”^) FR-II galaxy 3C 111 shows 
the differences between these two types of sources. A 
larger Doppler factor S = 14, a smaller emitting region 
(B = 2 X 10^® cm) and larger magnetic field (B = 0.04 G) 
are used to model 3C 111 ( de Jong et al.|2012| . The smaller 
volume of the emitting region used to model 3C 111 com¬ 
pared to M 87 yields a lower flux, but the strong Doppler 
factor S and larger magnetic field of 3C 111 increase the flux 
strongly, consistent with the larger bolometric luminosity of 
this source compared to M 87. 

As FR-I radio galaxies like M 87 are the parent popu¬ 
lation of BL Lacs, comparing the resulting SED of M 87 
with a BL Lac object Mrk 421 will illustrate the differ¬ 
ences between these two source classes. Using a one-zone 


SSC model, Blazejowski et al. (20051 found that Mrk 421 
(Lboi = 3 X lO^ergs”^) requires a beaming factor of 5 = 10, 
an emitting region a radius of B = 7.0 x 10^® cm and mag¬ 


netic field B = 0.405 G. The stronger beaming, due to the 
small jet angle of Mrk 421, boosts the intrinsic emission 
strongly. While the emitting region is smaller compared 
to that used to model M 87, the magnetic field is much 
stronger, in combination increasing the flux. 


4.4 M 87, a radio galaxy with a low luminosity 
BL Lac core 

In the overall context of SED models of y-ray bright sources, 
we can compare the derived values with the average ones 
for Fermi/LAT detected BL Lacs and FSRQs. As the FR-I 
radio galaxies can be understood as the parent population 
of the BL Lacs (e.g. Urry &: Padovani||1995 1, also their SED 
parameters are closer to the ones found for BL Lacs than 
for the brighter FSRQ class. As pointed out by |Ghisellini| 


et al. (20101, the BL Lacs on average appear to have similar 
masses as the FSRQs, with values around 10® — 10® Mq. 
The SMBH in M 87 exceeds this average value, and the 
bolometric luminosity Bboi — lO^'^ — 10'*® erg s”'^ results in 
an Eddington ratio of only A = Bboi/LEdd — 10“® — 10“® 
( Owen et al.||20(i0 Di Matteo et al. 20031. Nevertheless, a 
jet model similar to BL Lacs can be applied to the M 87 
SED. 

Only 15 radio galaxies and a few misaligned blazars 
have been detected with Fermi/LAT, compared to more 
than 1500 observed blazars ( Ackermann et al.|2015 (. In the 
case of blazars the y-ray emission is postulated to originate 
in the relativistic jet, which is pointed towards the observer. 
For M 87 the observed jet angle is 9 — 15°, causing the 
emission to be de-boosted as shown by the low Doppler fac¬ 
tor {S ~ 3.9 ) of the emission. Even though the origin of the 
y-ray emission is not completely clear, the modelling with a 
SSC model has shown that the high-energy emission is likely 
to arise in the jet. 

BL Lac sources are thought to have an advection domi¬ 
nated accretion flow (ADAF; [Reynolds et al.||1996| jPi Mat- 


teo et al. 20031 rather than an accretion disk, which is 


present in FSRQs. Due to the ADAF, the accretion is less 
efficient, and as such causes the source to be less luminous 
compared to FSRQs. Since in the unified scheme FR-I ra¬ 
dio galaxies, such as M 87 are considered to be misaligned 
BL Lacs, FR-I galaxies are thought to be powered by an 
ADAF as well. In the case of M 87 the broad-band emission 
is shown to be jet-dominated. 

The properties of M 87’s central engine are closer to 
that of BL Lacs than to FSRQs. The strength of the mag¬ 
netic field, the bulk Lorentz factor of the jet and also 
the electrons’ Lorentz factor distribution indicate rather 
a BL Lac type emission. Also the ov erall bolometric lu- 
minosity of Bboi — 10** — 10*®ergs“* 


(Owen et al. 


2000 


jPi Matteo et aT][^ giving an Eddington ratio of only 
A ~ 10~° —10~® points in this direction. The observed phys¬ 
ical properties of this AGN put it in a low luminous BL Lac 
class. The low power jet can be explained by the low ac¬ 
cretion rate, rather than by the large mass, although the 


finding by Laor (20001 that AGN with black hole masses 
larger than Mbh > 10® Mq are radio-loud, seems to hold 
also in the case of M 87. 

3C 120, another FR-I radio galaxy, was not included in 
the first, second or third Fermi catalogue, but using a 15- 


month data set Abdo et al. (2010c I derived a significance 
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of ~ 5.6cr between 0.1-100 GeV for this source. This source 
undergoes series of flares with a low long-term average flux 
(e.g. Sahakyan et al.||2015 |. 3C 120 has been monitored in 
the radio, UV and X-ray band ( Lohhnk et al.|[2013 K From 
the X-ray observations by Suzaku and XMM-Newton a 0.4- 
10 keV spectrum is derived, where a fit with a composite 
jet-baccretion disk model is favoured. The broad-band SED 
of this source has shown that the jet dominates the radio 
and 7 -ray emission, and contributes only ~ 10 % in the opti¬ 
cal, UV and X-ray bands (Kataoka et al.|[2011|. The 7 -ray 


to be a weak BL Lac object, consistent with the advection 


dominated accretion flow model for the core (Ptak et al. 


spectrum is soft, with F-, ~ 2.7, and the source is not de¬ 
tected in the TeV band, implying that the inverse Compton 
peak of this source is located in the X-ray to MeV band. 
In the case of 3C 120, as the 7 -ray emission is variable on 


the time scale of months. Sahakyan et al. (20151 argue that 
the emission comes from a compact, relativistically moving 
emission region, but exclude the jet knots as the origin as 
these are deemed too large. 

Only a few of the 7 -ray detected radio galaxies have 
been observed in VHE. Recently the MAGIC telescopes have 
detected another radio galaxy, the FR-I IC 310 at energies 
E > 300 GeV ( [Aleksic et al.|[^14[ ) . The SED of this source 
shows an inverse Compton peak in the TeV range. Com¬ 
bined with the low luminosity of this source, it was argued 
that IC 310 is an extreme case within the blazar sequence, 
with an extremely low accretion rate. Another possibility 
is that IC 310 is rather a misaligned version of an extreme 
BL Lac, where the low luminosity is connected to the large 
viewing angle. This does, however not explain the observed 
TeV variability. M 87 is also observed in the VHE band, and 
shows an inverse Compton peak in the hard X-ray/soft 7 - 
ray band. With its low bolometric luminosity, M87 is closer 
to IC 310 than to 3C 120, since both IC 310 and M87 have 
been detected firmly in both the 7 -ray and TeV band. 

The low Eddington ratio can be understood considering 
that the M 87 core is likely not fed by an accretion disk but 
by a radiatively inefficient accretion flow (RIAF) or advec¬ 
tion dominated accretion flow (ADAF; Reynolds et al.|1996 


Di Matteo et al.|2003 1. The assumption that the accretion is 


radiatively Inefficient also explains why in the case of M 87 
we do not see a significant thermal inverse Compton compo¬ 
nent in the X-rays, and an optical spectrum consistent with 
that of a LINER. This is different to other 7 -ray bright radio 
galaxies, which show an optical Seyfert core as e.g. Cen A 
( Beckmann et alT]|2011 l and 3C 111 ( de Jong et al.|2012 l. 


1998 Reynolds et al. 1996 \ and the overall low-luminous FR- 


I nature of this galaxy. The high X-ray flux detected with 
Suzaku at the end of 2006 seems to indicate the source was 
undergoing an outburst or flaring episode. The steep slope 
of the spectrum, with a power law index of F = 2 . 81 q '4 
tween 20-60 keV, indicates that the emission was likely the 
high-energy tail of the synchrotron branch. Using simulta¬ 
neous observations we created an SED for both the core and 
the jet knot HST-1 which is known to flare in other wave¬ 
bands such as the radio and optical. Due the non-imaging 
character of the Suzaku/PIN observations, we are not able 
to determine whether the enhanced emission results rather 
from the core of M 87 or in the jet based on the X-ray ob¬ 
servations alone, but from the SED modelling we conclude 
that the jet knot is the more likely candidate for the hard 
X-ray emission detected in 2006. 

In the unification model for AGN, radio galaxies are the 
counterparts of blazars, where the lower luminosity BL Lacs 
are linked to FR-I galaxies and the more powerful flat spec¬ 
trum radio quasars (FSRQ) to the FR-II sources. One way 
to account for the differences in luminosity and Compton 
dominance between the two source classes is to model the 
SED of BL Lacs with a simple SSC model and the FSRQ 
with a more complex model, for example using an external 


Compton component in addition to the SSC (Ghisellini et al. 
|1998p . The simple SSC model for the overall SED appears to 
be valid for this class of Fermi/LAT detected radio galaxies. 
Also in the cases of the brighter FR-II objects, like 3C 111, 
no external Compton component seems to be necessary to 
represent the SED, which is not in line with the unifica¬ 
tion model. Thus, in all these cases the dominating emission 
region is either far from a strong field of external seed pho¬ 
tons, like the broad line region (as in the case of 3C 111), 
or the broad line region itself is weak because of radiatively 
inefficient accretion, as might be the case in M 87. The hy¬ 
pothesis, that EC is not significant in 7 -ray detected radio 
galaxies should be tested in the case of the Fermi/LAT de¬ 
tected steep spectrum radio quasar 3C 207.0, which hosts a 
Seyfert 1.2 core. At a redshift oi z = 0.68 this object has a lu¬ 
minosity of L(2—10 keV) = 2.3x 10^® ergs“^, and the strong 
Seyfert core that displays an iron Ka line with EW ~ 60 eV 
should give rise to a significant photon field able to provide 
ample seed photons for inverse Compton processes in this 
case. 


5 CONCLUSION 

We report, for the first time, a hard X-ray detection 
of the FR-I radio galaxy M 87 using Suzaku/PIN data. 
The observations were made between November 29 and 
December 2, 2006 with an elapsed time of 187 ks, re¬ 
sulting in a flux of / = 10 “^^ erg cm“^ s“^ between 20 
and 60 keV. In addition, we derive a 3a upper limit of 
/(20 — 60keV) < 3 x 10“'^^ ergcrn”^ s“^ for the multi-year 
time averaged emission, based on 1.7 Ms of INTEGRAL 
IBIS/ISGRI data. 

By modelling the broad-band energy distribution with a 
one-zone SSC model we connect the average hard X-ray up¬ 
per limit to the emission and 7 -ray emission to the core emis¬ 
sion. The SED parameters show that M 87 can be considered 
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